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LONG-TERM  GOALS 

In  the  final  year  of  this  project,  the  primary  objective  was  the  application  of  the  propagation  model  to 
synthetic  aperture  sonar  simulations  and  simulations  for  scattering  within  a  homogeneous  waveguide. 
Refinements  to  image  source  enumeration  were  implemented  to  ensure  that  time-of-flight  wave 
packets,  associated  with  an  image  source,  arrive  at  a  receiver  location  within  a  time  window  of  interest. 

OBJECTIVES 

The  wave  propagation  model  developed  during  FY1 1-12  was  applied  to  several  scattering  problems, 
which  included  both  waveguide  simulations  (i.e.,  reflections  from  upper  and  lower  boundaries)  and 
reduced  half-space  simulations  (i.e.,  only  interactions  with  a  water-sediment  interface).  For  the 
scattering  problems,  there  are  two  wave  propagation  processes:  (1)  propagation  from  an  image  source 
to  a  target  location,  and  (2)  propagation  from  the  target  location  to  an  image  receiver.  During  step  (1), 
the  target  location  acts  has  a  virtual  receiver,  and  after  the  scattering  process  the  target  location 
becomes  a  virtual  source.  Thus,  the  refinement  to  the  image  source  enumeration  applies  to  the  image 
receiver  enumeration. 


APPROACH 


A  ray-base  model  developed  in  FY 1 1  assumed  a  homogeneous  layer  of  water  between  an  upper  semi¬ 
infinite  half-space  of  air  and  a  lower  semi-infinite  half-space  of  a  homogenous  sediment.  Rays  are 
assumed  to  travel  in  straight  line  segments  from  a  source  location  to  a  target  location  and  after 
interacting  with  the  target  to  a  receiver  location.  Within  the  current  report,  the  model  is  restricted  to 
the  propagation  from  an  image  source  to  a  target  location  (i.e.,  a  virtual  receiver).  For  an  interaction 
of  a  ray  with  an  interface,  an  appropriate  reflection  coefficient  is  included  in  the  received  signal.  At 
present  the  sediment  can  be  modeled  as  either  an  attenuating  fluid  with  a  frequency-independent  loss 
parameter  or  an  effective  density  fluid  model.  For  an  image  source  located  at  r,  and  a  target  located  at 
rt,  the  horizontal  and  total  separation  distances  are  | Rt  -  R\  and  \rt  -  r\,  respectively.  The  contribution 
of  the  z'th  source  to  the  frequency  spectrum  at  the  target  can  be  written  as 


=  PsrcW’  o 


exp  (iojtj) 


(1) 
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where  Psrc( co)  is  the  frequency  spectrum  of  the  transmitted  wave  packet  and  ro  is  a  reference  distance 
used  in  the  calibration  of  the  source  level.  The  quantity  within  the  square  brackets  accounts  for  the 
propagation  of  the  ray.  The  reflection  coefficients  at  the  upper  and  lower  boundaries  are  A(0,)  and 
B(Qj).  The  local  grazing  angle  0,  at  a  boundary  is  given  by  cos(0,)  =  \Rt  -  R, ]  /  | rt  -  r,|,  and  the  m  and  n 
exponents  indicate  the  number  of  interactions  a  ray  has  with  a  given  boundary.  The  time  delay  for 
propagation  from  the  source  to  the  target  is  tt  =  \rt—  r,\/c.  Here,  c  is  the  speed  of  sound  in  the 
waveguide.  If  the  scattering  from  the  target  to  the  /th  receiver  were  to  be  simulated,  then  Eq.  (1)  would 
be  multiplied  by  the  target’s  free-field  scattering  amplitude /(0,y,(j)//,co)  and  a  second  factor  similar  to 
the  quantity  in  the  square  brackets  with  t—*j  and  i—*t.  The  target-centered  polar  and  azimuthal  angles, 
Qij  and  (j);/  are  related  to  the  local  grazing  angles  0,  and  07. 

WORK  COMPLETED 

During  FY13,  simulations  of  scattering  from  a  15-m  radius  fluid  sphere  in  a  4000-m  deep  waveguide 
with  source/receiver/target  separation  distances  on  the  order  of  100  km  revealed  the  implementation  of 
Eq.  (1)  violated  reciprocity  under  certain  conditions.  If  the  source  and  receiver  were  located  below  the 
mid-depth  of  the  waveguide  the  implementation  satisfied  reciprocity.  If,  however,  either  the  source 
location  or  receiver  location  or  both  locations  were  above  the  mid-depth  of  the  waveguide  then 
reciprocity  failed. 

The  failure  was  caused  by  a  change  in  the  enumeration  of  image  sources,  the  method  used  to  detennine 
the  reflection  coefficient  exponents,  and  the  selection  of  image  sources  that  contribute  within  a  time 
window  of  interest.  The  depth  of  the  waveguide  is  h,  and  the  actual  source  is  located  at  zs  (>  0)  with 
the  water-sediment  boundary  at  z  =  0.  Previously,  the  first  image  source  in  a  set  of  image  sources  was 
placed  below  the  water-sediment  boundary  with  z\  =  -  zs,  and  the  c/th  quartet  of  image  sources  was 
added  to  the  set  of  image  sources  with  the  following  z-coordinates:  2 qh  +  zs,  2 qh  -  zs,  ~2qh  -  zs,  and 
-2 qh  +  zs  with  q  =  1 ,  2,  3,  ....  The  initial  ordering  assumed  that  zs  <  hi 2,  and  the  detennination  of  the 
reflection  coefficients  assumed  that  the  ray  associated  wiith  the  first  image  source  always  interacted 
with  the  lower  boundary  and  the  rays  in  the  enumeration  of  the  remaining  image  source  then  alternated 
with  interacting  with  the  upper  and  lower  boundaries,  respectively.  As  long  as  zs  <  hi 2,  z\  =  -  zs,  and 
the  addition  of  quartets  followed  the  prescribed  ordering,  then  the  image  sources  were  sorted  in  a 
manner  consistent  with  the  detennination  of  the  m  and  n  exponents  in  Eq.  (1),  and  reciprocity  was 
ensured. 

In  the  FY13  simulations  with  separation  distances  on  the  order  of  100  km,  it  was  advantageous  to 
consider  only  a  finite  length  time  window  with  a  defined  time  offset.  It  then  followed  that  only  image 
sources  with  propagation  delays  tt  =  | r,  -  r\lc  within  the  time  window  needed  to  be  considered;  and  so, 
the  set  of  image  sources  were  sorted  by  separation  distances.  When  zs  >  hi 2,  the  new  sorting  led  to  a 
first  image  source  with  z\  =  2 h  -  zs,  and  the  ray  associated  with  it  reflects  initially  from  the  air-water 
interface.  The  implementation  of  the  new  model  failed  to  updated  the  determination  of  m  and  n,  and 
thus,  reciprocity  was  violated. 

RESULTS 

As  an  example  to  demonstrate  that  reciprocity  holds  for  Eq.  (1),  a  wave  packet  with  a  10  Hz  carrier 
frequency,  6  ms  duration  and  Gaussian  envelope  was  propagated  from  a  source  to  a  receiver  in  an 
homogenous  water  waveguide.  The  water  depth  was  4000  m  with  a  sound  speed  and  density  of  c\  = 
1500  m/s  and  pi  =  1000  kg/m  .  The  sediment  properties  were  ci  =  1600  m/s,  p2  =  2000  kg/m  ,  and  8  = 
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0.003  with  8  being  the  loss  parameter.  The  source  and  receiver  were  initially  placed  at  (0,  0,  1)  and  (0, 
100,  1)  with  units  of  km.  As  a  comparison,  the  wave  packet  was  propgated  by  the  ray  model  of  Eq.  (1) 
and  a  model  based  on  the  trapped  modes  in  a  waveguide.  Figure  la  and  lb  show  the  case  where  the 
source  and  receiver  are  below  the  mid-depth  and  Figures  lc  and  Id  show  the  case  where  the  source 
(receiver)  is  below  the  mid-depth  and  the  receiver  (source)  is  above. 


Figure  1:  (a)  Propagation  for  a  source  (receiver)  and  a  receiver  (source)  below  the  mid-depth,  (b) 


Close-up  of  propagated  signal  in  (a)  demonstrating  reciprocity,  (c)  Propagation  for  a  source 
(receiver)  below  the  mid-depth  and  a  receiver  (source)  above  the  mid-depth,  (d)  Close-up  of 
propagated  signal  in  (c)  demonstrating  reciprocity.The  location  of  the  points  are  A  =  (0,0,1)  km,  B  = 
(0,100,1)  km,  C  =  (0,  0, 1)  km,  and  D  =  (0, 100,  3)  km.  Note,  the  red  and  green  curves  overlap  and 

the  blue  and  green  curves  overlap. 


The  discrepancy  between  the  ray  model  and  trapped  mode  model  has  been  traced  to  an  unphysical 
discontinuity  in  the  frequency  of  the  Green’s  function  generated  with  the  trapped  modes  model.  The 
cause  of  this  discontinuity  is  currently  unknown. 
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RELATED  PROJECTS 


The  companion  award  (N00014-07-G-0557,  Task  0023),  “Generation  of  Synthetic  SAS  Data  for 
Targets  near  the  Seafloor:  Target  Scattering  Component,”  seeks  to  improve  scattering  models  for 
targets  near  an  interface  used  in  SAS-ST.  Improvements  to  the  propagation  model  enhance  the 
combined  propagation  and  target  scattering  model. 

The  results  of  this  project  has  been  combined  with  the  target  scattering  component  and  used  as 
benchmark  for  the  DARPA  project,  “Exploiting  Seismic  Events  to  Maintain  Under-Ice  Awareness,” 
(Award  Number  W91 1QX-12-C-0163).  The  project  has  developed  a  scattering  model,  based  on  the 
trapped  modes  in  an  inhomogeneous  waveguide. 

The  research  of  this  project  in  combination  with  the  target  scattering  component  has  been  used 
extensively  in  data-model  comparison  for  the  ONR  funded  project,  “Acoustic  Color  of  Mines  and 
Mine-like  Objects,”  (N00014-07-G-0557/0032).  It  provides  a  fast  robust  model  for  the  simulation  of 
target  scattering  near  and  at  an  interface. 
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